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Abstract 

A  kinetic  study  is  made  of  a  system  of  neutral 
particles  which  interact  under  gravitational  forces. 
The  Jeans  instability  is  recaptured  in  the  limit  of 
long  wavelengths  where  the  growth  rate  assumes  a 
characteristic  maximum  value »   As  the  wavelength 
diminishes  to  a  certain  critical  distance,  the  insta- 
bility vanishes.   This  distance  is  consistent  with 
thermal-gravitational  estimates  for  the  diameter  of  a 
star.   Plasma  oscillations  are  found  in  the  central 
region  of  a  one-component  charged  gas  cloud. 
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Gravitational  Instability  and  One -Component  Plasma  Oscillations 

I .   Introduction  and  Summary  of  Results 

The  equations  which  govern  a  "colllsionless"  one-component 
gas  which  is  composed  of  neutral  particles  are  very  similar  to 
those  which  govern  a  "colllsionless"  plasma.   The  single  dif- 
ference resides  in  the  change  of  sign  between  the  two  divergence 
equations ,    viz , , 

(1)  Natural    V'S,  =  -  ctp 

(2)  Charged    V-E^  -  +  Pq   - 

The  gravitational  force  field  is  G^,  and  E__  is  the  electric 
field.  The  mass  density  is  p,  q  is  charge  density,  and  a  and 
P   are  positive  constants. 

In  this  paper  the  consequences  of  this  single  difference 
are  examined.   Although  the  plasma  problem  has  been  extensively 
studied  in  the  past  score  of  years  (Landau  19-^6,  Van  Kampen  1957)  ? 
interest  in  the  neutral  gas  problem  has  only  recently  been  revived 
(Lynden-Bell  I962,  R.  Simon  1962 ) ,  having  first  been  examined  in 
the  macroscopic  limit  by  Jeans  in  1902=   The  present  work  follows 
quite  closely  that  of  Simon. 

The  first  observation  is  that,  while  the  roots  of  the  Landau 
problem  are  obtained  only  by  affecting  an  analytic  continuation, 
the  relevant  solutions  of  the  gravity  problem  are  strictly  normal 
modes.   Secondly,  the  neutral  gas  modes  exhibit  pure  instability 
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as  opposed  to  the  fact  that  all  solutions  of  the  plasma  problem 
are  damped.  In  addition,  all  of  the  gravity  modes  which  propa- 
gate are  damped. 

The  growing  fluctuations  are  examined  in  detail.   The 
maximum  growth  rate  of  these  instabilities  is  found  to  be 

CO  =  l/amn   ,  where   m   is  the  mass  of  the  elements  and   n   the 
O    '    o  '  o 

equilibrium  number  density.   This  maximum  growth  rate  is  asso- 
ciated with  very  long  wavelength  disturbances  --  which  is  in 
agreement  with  Jeans'  (1902)  conclusions.   The  growth  rate 
decreases  to  zero  as  the  wavelength  of  the  disturbance  approaches 
the  distance   d  =  / C /cu    (Simon  I962).   The  mean  thermal  velocity 
of  the  equilibrium  distribution  is   C.   The  distance   d   is  the 
analogue  of  the  Debye  (1923)  distance  for  plasmas.   It  is  also 
the  well-known  lowest  order  estimate  for  the  dimensions  of  a 
star  --  which,  in  turn,  is  obtained  by  equating  the  random  thermal 
energy  of  a  surface  particle  to  its  gravitational  potential 
energy. 

Finally,  a  one-component  charged  spherical  plasma  is 
examined,  and  it  is  found  that  plasma  oscillations  exist  in  the 
central  region.   Such  fluctuations  should  exhibit  themselves 
in  a  non-charge-neutral  universe  (Lyttleton  and  Bondi  1959?  and 
Hoyle  i960). 
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II .   Analysis 


Let  us  consider  a  neutral  gas  which  is  in  a  Maxwellian 
state   f  ,  which  exists  in  the  absence  of  any  inter-particle 
interaction  for  times  t  <  0.   At  t  =  0,  the  interaction  G 
is  "turned  on".   The  subsequent  distribution  is   f  +  g,  where 
G   and   g   are  first  order  infinitesimals.   Neglecting  second 
order  terms  yields  the  following  equations  for  g   and  G: 


O) 


+ 


^V^ 


e^-^^^o 


The  linearization  of  (l)  gives 


(^) 


U*G  =  -  amn   -  am 
V.  _        o 


gd  ^ 


The  mass  of  the  particles  is   m,  and  the  distribution  function 
f    is  normalized  over  velocity  I   to  the  number  density  n  . 


If  the  plane  wave  forms,   exp  <  -icot  +  ik'2S.  '    ^^^    substi- 
tuted into  the  above  set,  there  results 


(5) 


am 


r  v.f  d^^ 

5l  o 


+   k 


'G 


lamn 


For  longitudinal  (k  X  G  =  O)  fluctuations,  the  following 
dispersion  relation  is  obtained  from  the  above: 


(6) 


r 


k 


V  (p  dv 
V  -  (co/k) 
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This  equation  Is  In  purely  non-dimensional  form.   The  (non-  ■ 
dimensional)  wave  number  k   is  written  for  kd,  the  frequency 
00  is  written  for  (oj/cd  ),  and  the  velocity  v   is  written  for 
i/C.      The  function   ())   is  the  one-dimensional  distribution 
function, 


(7)     ^(v)  =  (27r)-^/2  e-^  /2  ^  (2^)-V2  ^-3 


d^  d^   e"^  /^^ 
X  ^y 


The  distance   d   is  the  analogue  of  the  Debye  distance  in  plasmas 
(p.  Debye  and  E.  HiUckle  I923),  while  the  growth  rate   co   is  the 
analogue  of  the  characteristic  frequency  of  plasmas  (L.  Tonks 
and  I.  Langmuir  I929): 


(8) 


00  =  amn 
o      c 


(9) 


^2    „2/  2 
d   =  C  /cu 

'  o 


Equation  (6)  will  have  roots  only  when  the  function 


(10) 


P(z)  = 


r   ({)dv 

v-z 


co/k 


Is  real  and  positive. 

It  is  easily  shown  (see  Appendix)  that   Im  P  >  0   in  the 
first  quadrant  of  the  z-plane  (including  the  positive  real 
axis),   Im  F  <  0   in  the  second  quadrant  (including  the  negative 
real  axis),  and  Im  F  =  0  along  the  imaginary  axis.   It  follows 
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that  all  the  unstable  roots  of  (6)  (i.e.,  roots  lying  In  the 
upper  half  plane)  lie  along  the  imaginary  axis,  or  equivalently 
for 


11 


z  =  IP  , 


For  these  values  of   z,   F  assumes  the  form, 

CO 

f        2 


(12) 


F(1P)  = 


J 


V  (p  gy 

V  +  P 


If   k   is  rewritten  as 


(13) 


,2     2/^2 
k   =  V  /p   , 


where   v   is  the  growth  rate. 


(14) 


CD  =  IV 


(all  values  are  non-dimensional  through   d   and  cu  ),  then  (6) 
now  appears  as 


(15) 


2/„2 
V  /p 


/    2  L 

V  (p  dv 

2   q2 

V  +  6 


For  a  given  frequency   v,  this  equation  determines   p   which, 
in  turn,  through  (13) j  gives   k.   Both  sides  of  (15)  are  plotted 
in  Fig.  1 . 


V>(jJ, 


F 


Figure  1 

The  Inequalities  which  are  depicted  In  Fig.  1  are  obtained 
by  examining  the  asymptotic  behavior  of  F  for  large  p, 
according  to  which 


(16) 
(13) 


F.-^  (1  -  ^  +  ...)  , 


k  =  V  /p   . 


It  follows  that,  for  large  P,  both  k   and  the  first 

2 
derivative  of  k   are  larger  than  the  related  values  of  F  if 


v"  >  1  , 
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which  appears  as 

(17)  V  >  CO 


o 


in  dimensional  form.   It  follows  that  there  are  no  solutions 
for   V  >  o)  ,  so  that   uj   represents  the  maximum  growth  rate  of 
these  Instabilities. 

The  relevant  root  for  large  p   is  obtained  by  equating 
equations  (15)  and  (l6).   There  results  (in  dimensional  form). 


(18)  v^  =  CD^  -  3k^C^ 

\    /  o 


which  is  identical  to  Jeans'  result.   These  values  of   v   are 
close  to  the  largest  (i.e.,   co  )  and  clearly  relate  to   k  ^  0, 
or,  equivalently,  to   A  ^=s^  oo  ,   which  is  again  consistent  with 
Jeans'  conclusions. 

The  smallest   v   values  lie  at  the  opposite  extreme, 
P  '^  0.   The  expansion  of  F(iP)   about  small   p   appears  as 


(19)         F  '^  1  -  p^p 


r 

dv  +  ...  -  p  iirpl  e 


i  dv  +  ...  -  rV^TTq   pP  /2 


2 

v 


=  1  -  ^Ytt/^  + 


The  symbol   P  :  indicates  that  the  principal  value  is  to  be 

2 
taken.   Equating  the  above  series  to   k   gives  the  desired  roots. 


(20)  k^  =  v^/p^  =1  -  p  +  ... 
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It  follows  that  in  the  limit   p  — >  0,   v  — >  0  and  k  — >  1 .   In 
dimensional  form. 


;2i)  V  «  cD^  ,    k  ->  d  ■•■ 


We  conclude  that,  although  long  wavelength  disturbances 
grow  most  rapidly,  this  growth  rate  diminishes  to  zero  as  the 
distance  between  condensations  approaches  the  critical  distance 
d.   In  a  homogeneous  gas  of  infinite  extent  the  characteristic 
time  for  the  development  of  this  final  state  is  cd 

All  of  the  remaining  modes  are  damped.   The  exact  forms  are 
obtained  by  continuing  F  analytically  into  the  lower  half 
z-plane.   The  only  physically  relevant  fact  about  these  damped 
modes  is  that  some  of  them  propagate  (viz.,  the  roots  that  lie 
off  the  imaginary  axis). 

In  the  one-component  charged  gas  problem,  the  E_  field, 
due  to  the  Coulomb  interaction,  is  very  much  larger  than  the 
Newtonian  interaction  force  encountered  above,  so  that  a  pertur- 
bation technique  such  as  was  used  in  the  neutral  gas  problem 
would  appear  inconsistent  in  the  charged  gas  case.   However,  as 
is  well  known,  the  Coulomb  force  becomes  vanishingly  small  in 
the  center  of  a  charged  spherical  cloud,  being  of  the  order  (r/R), 
where   R   is  the  radius  of  the  cloud  and   r  the  radius  of  the 


* 

2  2 

That   k   lies  in  the  range   0  <  k  <  1   is  most  easily  seen 

by  plotting  k  -   constant  vs.  P  (cf.  Fig.  l). 
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said  central  region.   The  related  equations  are  identical  to 
the  set  (2,  3)  with   G^  replaced  by  E_^  and  a   by   -  p.   As 
Landau  (19^6)  has  shown,  all  solutions  are  damped.   In  the 
limit  of  long  wavelengths,  the  damping  is  exponentially  small, 
and  the  roots  appear  as 


(22)  CD  -  cDp  +  3k  C   , 


a  propagating  oscillatory  mode.   Although  these  modes  are  quite 
well  known,  they  are  usually  obtained  assuming  charge-neutrality 
in  the  equilibrium  configuration.   The  classical  plasma  fre- 
quency is  cOp  (Tonks  and  Langmulr  1929)' 
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Appendix 

The  following  arguments  exhibiting  the  properties  of  P(z) 
are  due  to  Professor  H.  Weltzner. 

First  we  recall  that   F(z)   is  defined  as 


(A.l 


F(z)  - 


V  e   ^   dv 
v-z 


If  the  Integrand  Is  multiplied  by   1  =  (v+z)/(v+z),  there  results 


(A. 2) 


F(z)  = 


2  -V 

V  e 

2     2  ' 

V  -  z 


where  the  remaining  term  is  dropped  because  of  oddness  of  the 

2   -2     2   -2 
Integrand.   Multiplying  again  by   1  =  (v  -  z  )/ (v  -  z  )   gives 

the  form 


(A. 3) 


P(z) 


V   e 


(v  -  z  ) 


v-z 


It  follows  that 


(A. 4) 


_2 
Im  F(z )  =  -  a  Im  z   , 


where  a   is  the  real  positive  number. 


(A. 5) 


a  ^ 


^   2  -v'^/2    , 
V   e   ^   dv 


V 


^^2? 


1  0 
Setting  z  =  re   ,  one  obtains 
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:a.6) 


Im  F  =  sin  2( 


ar 


For  values  of   z   along  the  real  axis,  F  appears  as 


(A. 7) 


F(x)  =  P 


'     -V  /2  ^  2/„ 

X£ :__^  +  27rix  e-^  /2 

v-x 


The  principal  part   P:   is  purely  real  while  the  Cauchy  contri- 
bution is  purely  imaginary,  so  that  for   9  =  0,  ir , 


(A. 8) 


Im  F  =  27r  X  e 


-X 


The  value  of   F   along  the  positive  imaginary  axis  has  been 
calculated  in  the  text.   It  follows  that  in  the  entire  upper 
half  z-plane. 


(A. 9) 


Im  F  >  0 
Im  F  =  0 
Im  F  <    0 


0   <    e   <    7r/2 
e   =   ir/2 
Tr/2   <   9   <   TT 


Q.E.D, 


-  14    - 


References 

1.  Debye,  P.  and  E.  Hi^ckle,  Physlk.  Z.  24_,  185  (1923). 

2.  Hoyle,  P.,  Proc .  Roy.  Soc .  A,  257,  ^31  (I96O). 
5.   Jeans,  J.,  Phil.  Trans.  I99,  A,  ^9  (1902). 

k.      Landau,  L. ,  J.  Phys.  (USSR)  10,  25  (19^6). 

5.  Lynden-Bell,  D. ,  Monthly  Notices,  Roy.  Ast .  Soc.  124,  279 
(1962). 

6.  Lyttleton,  R.  A.  and  H.  Bond!,  Proc.  Roy.  Soc.  A,  252,  313 

(1959). 

7.  Simon,  R.,  Bull,  de  I'Acad.  Royal  de  Belglque,  Series  5? 
47,  7  (1962). 

8.  Tonks,  L.  and  I.  Langmuir,  Phys.  Rev.  53,  195  (1929). 

9.  Van  Kampen,  N.  G.,  Physica  XXIII,  64l  (1957). 


-  15  - 


f  EB  1 3  1963 

DATE  DUE 


1 

/ 

/ 

CAYLO^D 

en    -.TEO                / 

/  nYu 

NYO- 

'   '^5|  -"ft; 

NYO~ 


c.l 


N.  Y.  U.  Courant  Institute  of 
Mathematical  Sciences 

4  Washington  Place 
New  York  3,  N.  Y. 


f 


